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ABSTRACT 


The  exact  inverse  scattering  integral  eqiation  of  this  author  is  applied  to  passive  seimic 
artillery  location  in  mlewan  media.  It  is  show  that  for  such  media  this  integral  equation 
need  not  be  solved,  since  a  tine-space  display  of  the  effectal  field  spatially  and  terporally 
locates  the  inpulsive  sources.  Niner ico-exper inenta)  verification  is  presented. 
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SECTION  I 


INTRODUCTION 


In  section  II  the  Rayleigh  surface  wave  equation  is  reformul ated  into  an  inhomogeneous 
free-space  wave  equation  with  a  source  term  depending  on  an  arbitrarily  chosen  reference 
velocity.  In  section  III  the  exact  inverse  scattering  integral  equation  of  this  author 
is  reder  ived  and  appl  i  ed  to  this  inhomogeneous  free-space  Rayleigh  surface  wave  equation. 
The  physical  meaning  of  the  effectal  field  is  discussed  in  section  IV.  Specifically,  that 
this  effectal  fi eld  consists  of  the  destructive  interference  between  the  true  causal  field 
and  a  fictitious  anti-causal  field,  both  propagating  at  the  reference  velocity,  and  that 
if  the  true  field  contains  hrpulsive  sources,  then  the  effectal  field  also  contains  inpuls  ive 
sources  at  the  same  spatial  and  temporal  locations.  A  graphic  time-space  display  of  the 
effectal  field  thus  yields  the  correct  temporal  and  spatial  locations  of  the  impulsive 
sources,  obviating  the  need  for  a  solution  of  the  inverse  scattering  integral  equation  for 
the  case  where  only  these  locations  are  needed,  and  the  unknown  medium  velocity  is  not 
needed.  In  section  V  the  two  dimensional  inverse  scattering  integral  equation  is  reformulated 
in  the  time  domain;  further  illustrating  the  physical  meaning  of  the  effectal  field.  Numerico- 
experimental  results  are  presented  in  section  VI. 

A  closed  form  solution  of  a  more  general  exact  inverse  scattering  integral  equation,  developed 
under  a  separate  contract,  is  presented  in  the  appendix  section  IX.  Numsr ico-exper  imental 
verification  of  this  solution  is  also  presented  in  section  VI. 

In  light  of  these  numer  ico-exper  imental  results,  reccrrmendat  ions  for  future  research  are 
made  in  sect  ion  VI I. 
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SECTION  II 


REFORMULATION  OF  THE  RAYLEIGH  SURFACE  WAVE  EQUATION 


The  two-dimensional  Rayleigh  surface  wave  equation  in  the  frequency  domain  is 


V2<ji+k^4i  =  0 


(1) 


where 


k 


v 
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v  (X,od) 
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(2) 


and  where  the  ve  I  oc  i  ty  j  $  related  to  the  shear  modulus  V*  and  the  mass  density0  by 


v 
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0.9194 


(3) 


Next,  let  a  potential  V  and  a  source  density  P  be  defined  respectively  in  terms  of 
an  arbitrarily  chosen  reference  velocity  c  as 


V(X,u>) 


(4) 


2 


p  (  X  ,  U)  ) 


V  (  X  ,  a) )  <f>  (  X  ,  w  ) 


(5) 


The  wave  equation  (1)  thus  becomes  the  inhomogeneous  wave  equation 


V  2  <f>  +  k2  <j>  =  -p 


subject  to  the  constitutive  equation  (4). 


(6) 
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SECTION  I  I  I 


THE  INVERSE  SCATTERING  INTEGRAL  EQUATION 


In  this  section  the  inverse  scattering  integral  equation  of  this  author  1 1  )  will  he 
rederived.  This  integral  equation  has  been  discussed  and  studied  extensively  by 
this  author  [2]  -  [7]  and  others  [8]  -  [28]. 

Let  a  field  *•  satisfy  the  inhomogeneous  wave  equation 

V2<p  +  k 2  #  =  -p  .  (7) 

Furthermore,  let  G  be  the  free-space  (relative  to  v  =  c )  Green's  function  for  (7), 
which  satisfies  the  inhomogeneous  wave  equation 

V2G+k2G=-6  »  (g) 

and  the  Somme rfeld  radiation  condition  at  infinity. 

next,  let  H  be  the  imaginary  part  of  the  Green's  function  G;  i.e., 

H  :  Im  G  (9) 


4 


It  thus  follows  from  (8)  that  this  imaginary  part  of  the  Green's  fund  i  on  sat  i  sf  ies 
the  homogeneous  wave  equation 


V2H  +  k2H  =  0 


(10) 


Next,  let  the  effectal  field  0  (the  term  "effectal"  will  be  explained  and  justified 
in  section  IV)  be  defined  as 


e 


dS  •  (  $>  VH 


H  74  ) 


(11) 


where  S  is  the  surface  on  which  the  field  4>  is  known,  and  where  the  support  of  the 
unknown  sources  P  is  inside  this  surface  S. 

By  Green's  theorem,  (11)  reduces  to 
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d  v  (  4>  V2H 


H  V2<p  ) 


(12) 


which,  by  (8)  and  (10),  further  reduces  to 
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(13) 


(14) 


The  unknown  sources  p  in  the  volume  V  are  thus  related  to  the  known  effectal  field  0  in 
the  volume  V  by  the  proper  (i.e.,  x  and  x'*V)  Fredholm  integral  equation  of  the  first 
kind;  where  the  effectal  field  0  in  the  volume  V  is  computable  from  knowledge  of  the 
field  4  on  the  surface  S  by  (11). 
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SECTION  IV 


THE  PHYSICAL  MEANING  OF  THE  EFFECTAL  F  I  E  L  n 


The  imaginary  part  of  the  Green's  function  can  he  written  as 


(151 


Since  the  Green's  function  is  causal,  it  follows  from  the  Fourier  transform  relationship 


G  (  r  ,  a;  r  /  c  )  q  (  r  ,  t  -  r  /  c  ) 


that 


* 

G  (  r  ,  ter  /n)  ci  (  r  ,  t  +r  /  e  1  .  (17) 


It  thus  follows  from  (15)  that 


H ( r,u' 
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In  the  time  domain,  the  effectal  field  ■'  is  thus  by  (14)  and  (18) 
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The  first  convolution  on  the  right  hand  side  of  (Id)  is  clearly  the  real  physical  t  ine- 
retarded  causal  field  0(x,t)  radiated  by  the  sources  hi  ( x ,  t  ) .  The  second  integral 
on  the  right  hand  side  of  (Id)  is,  however,  a  fictitious  time-advanced  anti-causal 
field  "radiated"  by  the  sources  n  (x  ,t  ).  The  effectal  field  "(x,t)  is  thus  an  imaeinarv 
field  due  to  the  destructive  interference  between  the  time-retarded  causal  field  and 
the  fictitious  time-advanced  anti-causal  field. 


Examination  of  (4),  (5),  and  (19)  readily  reveals  that  if  the  field  contains  spatially 
and  temporally  impulsive  sources,  then  the  effectal  field  will  also  contain  these  impulsive 
sources  at  the  correct  spatial  and  temporal  locations.  A  graphic  time-space  display  of 
the  effectal  field  will  thus  yield  these  correct  spatial  and  temporal  locations,  obviating 
the  need  for  a  solution  of  the  inverse  scattering  integral  equation  (19)  for  the  case  where 
only  these  locations  are  needed. 


SECTION  V 


THE  TWO-DIMENSIONAL  EFFECTAL  FIELD  IN  THE  TIME  DOMAIN 


In  the  frequency  domain,  the  two-dimensional  Green's  function  is 
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4  [  ° 


(kr)  +  iY  (kr) 
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(20) 

(21  ) 

(22) 


Thus,  in  the  frequency  domain,  the  imaginary  part  of  the  two-dimensional  Green's  function  is 


H  -  Tm  C- 


(23) 
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In  the  time  donain,  the  "imaginary  part"  of  the  two-dimensional  Green's  func t ion  becones 
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which,  with  the  aid  of  the  integral  representation  of  the  Vessel  functions 
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yields 
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(27) 


which,  in  turn,  with  the  aid  of  the  delta  function  representation 


6(t) 


(28) 


reduces  to 
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Introducing  the  change  of  variable 


z 


r 

—  cos 
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(30) 


yields  for  (29) 


h 


(31) 


whi  ch  yields 
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(32) 


It  is  now  of  interest  to  compare  (32)  with  the  full  two-dimensional  Green's  function  in 
the  time  domain  given  by  | 29  1 

0 

In  the  time  donain,  the  two-dimensional  inverse  scattering  integral  equation  (14)  is  thus 

n(X-X'.t-t')  p(X',t')  dt'  dv'  =-  e(X,t)  .  (34) 

As  previously  indicated,  examination  of  (4),  (5),  and  (34)  readily  reveals  that  if  the 
field  contains  spatially  and  temporally  impulsive  sources,  then  the  effectal  field  will 
also  contain  these  impulsive  sources  at  the  correct  spatial  and  temporal  locations. 
A  graphic  time-space  display  of  the  effectal  field  will  thus  yield  these  correct  spatial 
and  temporal  locations,  obviating  the  need  for  a  solution  of  the  inverse  scattering 
integral  equation  (34)  for  the  case  where  only  these  locations  are  needed. 
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SECTION  VI 


NUMER  ICO-EXPER  (MENTAL  RESULTS 


A  two-dimensional  surface  with  a  unity  normalized  Rayleigh  surface  wave  velocity  was 
selected.  A  spatially  and  temporally  impulsive  source  was  placed  in  this  surface, 
and  an  analytic  expression  for  the  radiated  field  was  obtained  for  all  points  on  the 
boundary  of  a  128x128  data  grid  containing  this  source.  The  effectal  field  at  all 
points  in  the  grid  was  numer i ca I  I y  computed  by  (11)  for  256  f r equenc i es ,  with  a  reference 
velocity  c=2  (i.e.,  100%  off  the  correct  velocity).  Next,  a  fast  Fourier  transform 
of  the  effectal  field  into  the  time  domain  was  taken.  A  graphic  time-space  representation 
of  every  fifth  point  for  t=-5  to  t=70  of  this  effectal  field  is  shown  on  page  13. 
As  shown  in  section  V,  the  correct  temporal  and  spatial  locations  of  the  impulsive  source 
are  obtained  (see  t=0  frame);  however,  as  anticipated,  the  effectal  field  "flies  in" 
and  "flies  out"  at  twice  the  correct  velocity. 

The  preceeding  numerical  experiment  was  repeated  with  a  random  zero  db.  signal  to  noise 
added  to  the  radiated  field  prior  to  the  computation  of  the  effectal  field  by  (11). 
A  similar  graphic  time-space  representation  for  this  computation  is  shown  on  page  13. 
Again,  the  correct  temporal  and  spatial  location  of  the  impulsive  source  is  obtained  (see 
t=0  frame),  and  again,  as  anticipated,  the  noise  corrupted  but  recognizable  effectal 
field  "flies  in"  and  "flies  out"  at  approximately  twice  the  correct  velocity. 

'Mien  this  research  was  initiated,  only  the  1973  inverse  scattering  integral  equation  of 
this  author  was  in  existence,  and  only  the  implementation  of  the  graphic  time-space 
representation  of  the  effectal  field  was  anticipated.  A  more  complete  inverse  scattering 
integral  equation  and  its  analytic  closed  form  solution  (presented  in  the  appendix) 
were  obtained  (under  separate  contract  N00014-76-C-0082 )  in  the  latter  part  of  this 
research.  This  solution  was  numer ico-exper imental ly  verified  for  a  point  source  and 


1 1 


a  point  source  pair  separated  by  1/5  of  a  wavelength,  under  the  same  conditions  as  the 
t  ime  -space  representations  of  pages  13  and  14,  except  that  the  implied  Hilbert  transform 
was  executed  with  the  aid  of  the  Wi  ener-Lee  transform  [30],  4  graphic  representation 
of  the  resulting  fields  is  on  pages  15  and  16. 
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SECTION  VII 


RECOMMENDATIONS  FOR  FUTURE  RESEARCH 


VWien  this  research  was  initiated,  only  the  1973  inverse  scattering  integral  equation  of 
this  author  was  in  existence,  and  only  the  implementation  of  the  graphic  time-space 
representation  of  the  effecta!  field  was  anticipated.  The  more  complete  inverse  scattering 
integral  equation  and  its  analytic  closed  form  solution  presented  in  the  appendix  were 
obtained  (under  separate  contract  N0001 4-76-C-0082 )  only  in  the  latter  part  of  this 
research.  The  successful  numer ico-exper imental  verification  of  this  new  closed  form 
solution  indicates  that  future  research  should  be  concentrated  toward  this  solution; 
specifically,  but  not  nescessarily  limited  to: 


1.  Investigate  and  determine  which  of  the  various  existing  numerical  Hilbert  transform 
and  Wiener-Lee  transform  algorithms  is  best  suited  for  the  solution. 

2.  Investigate  and  determine  the  effects  of  reference  slowness  sampling  density  and 
doma i n  t  runcat i on. 

3.  Investigate  and  determine  the  effects  of  incomplete  aspect  angle  information  and 
and  input  data  sampling  density. 

4.  Investigate  and  determine  the  effects  of  input  data  noise  and  errors. 

5.  Make  use  of  the  availability  of  wide  frequency  band  and/or  time  domain  information, 
for  the  purpose  of  overdetermining  the  problem  and  reducing  the  effects  of  items  1  through  4. 

6.  Make  use  of  item  5  when  a  priori  information  exists  that  the  medium  is  non-di  spers  ive. 
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EXACT  INVERSE  SCATTERING  THEORY 
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this  author's  report  "Exact  Inverse  Scattering  Theory",  March  1980,  vrtiich 
was  prepared  for  and  supported  by  the  Office  of  Naval  Research,  800  North 
Quincy  Street,  Arlington,  Virginia  22217,  Under  contract  NOOOI 4-76-C-0082. 
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APPENDIX  ABSTRACT 


The  concepts  of  reference  wave  slovmess  (reciprocal  of  velocity)  and  an 
associated  free  reference  space  Green's  fuction  slowness  spectrum  are 
introduced.  A  modified  Kirchhoff  surface  integral,  containing  only 
the  imaginary  part  of  this  free  reference  space  Green's  function  slowness 
spectrum,  is  formulated,  yielding  an  integral  equation  for  the  unknown 
fields  and  sourcei  in  the  interior  of  a  closed  surface  on  which  the 
(remotely  sensed)  fields  are  known.  Awell-oosed,  analytic  closed  form 
solution  of  this  integral  equation  is  obtained. 


I  NT  RODUCT I  ON 


Presented  is  a  unified  approach  and  solution  to  the  inverse  scattering  and  inverse 
source  problems  for  the  inhomogeneous  scalar  wave  equation 


ri2  x  ,  W 
V<  +  ? 


subject  to  the  constitutive  equation 


p  =  V  4>  , 


and  the  homogeneous  scalar  wave  equation 


V2  tp  +  -  <j>  =  0 

C2(X,w) 


To  this  end,  the  single  mixed  scalar  wave  equation 


V2  <p  +  -  4>  =  -  p 

c2(x,w) 


(2) 


m 


is  introduced.  From  an  inverse  scattering  inverse  source  perspective,  (4)  reduces 
to  (1 )  if  the  medi um  wave  velocity  c(X,w)  is  a  knovwi  constant  and  the  source  •’  is  the 


unknown,  and  (4)  reduces  to  (3)  if  the  sources  are  known  to  he  zero  and  tie  medium 
wave  velocity  :x,.)  is  the  unknown. 

It  is  argued  that  the  inverse  solution  presented  is  an  alternative  f  to  the  lirei.? 
Kirchhoff)  integration  of  the  wave  equation.  It  is  thus  appropriate  to  t  -vi  -., 
relevant  properties  of  the  direct  Kirchhoff  integration 

:  j  .  fs- 

of  the  wave  equation  (1  ).  Specifically,  the  surface  integral  in  (SI  is  to  e.iu  v  t  !  v  n 
st  it  arent  relating  the  field  at  a  field  paint  on  one  s  i  ic  of  the  closet  surf  ur  nr  ■  i  .  ‘ 
by  all  the  sources  on  the  other  side  of  the  closed  surface,  via  the  fields  produced  ‘v 
these  sources  on  this  closed  surface.  The  inverse  scattering  inverse  source  problem 
is,  however,  characterized  by  both  the  field  point  for  the  unknown  fields  as  veil  as  i1 ! 
the  unknown  sources  being  on  the  same  side  of  the  closed  surface  (on  which  the  r*";ir 
sensing  is  accompl  i shed) ,  for  which  situation  the  Kirchhoff  surface  integral  vanishes, 
thus  rendering  this  Kirchhoff  surface  integral  useless  for  the  inverse  scattering 
inverse  source  probl cm.  A  niodified  Kirchhoff  surface  integral,  which  does  not  suffer 
from  this  pathology,  is  introduced  next. 


2  S 


SECT  ION  A  - I  I 


THE  INVERSE  SCATTERING  INTEGRAL  EQUATION 


Let  G  be  the  free  reference  space  Green's  function  satisfying  the  inhomogeneous 
wave  equation 


(«> 


and  the  Sommerfeld  radiation  condition  at  infinity,  where  v  is  any  arbitrarily  chosen 
reference  velocity. 

Next,  let  an  effectal  field  0  be  defined  as 


5 


(71 


where 


G .  =  Tm  G  , 

which,  by  (6),  satisfies  the  homogeneous  wave  equation 


VZG.  +  0.0 

t  v'  i 


(8) 


(9) 
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By  Green's 


which,  by 


theorem,  (7)  reduces  to 


J 


6  =  I  dv  ($  v2g.  -  G^  V2$) 


(4)  and  (9),  reduces  to 


0  -  ^  dv  C<J)(~~2  G_. )  -  Gv  ( - tp~o  )  ]  , 


V  *  1  C2  ( X  ,  u) ) 


r  r  r  u2  11)2 1 

J  J  I  c2 (X,ui)  v2  J 

w  \/  ** 


(10) 


(11) 


(12) 


(which,  for  the  case  of  a  known  constant  medi  im  wave  velocity  c,  unknown  sources;  , 
and  a  reference  velocity  chosen  as  v=c ,  reduces  to  the  earlier  inverse  scattering 
integral  equation  of  this  author  [1]). 


SECT  ION  A- I  I  I 


SOLUTION  OF  THE  INTEGRAL  EQUATION 


Let  a  medium  reference  wave  slowness  ^  be  introduced  and  defined  as 


o 


V 


(13) 


The  inverse  scattering  integral  equation  (12),  in  terms  of  this  reference  slowness, 
thus  is 


e(x,u),o) 


V 


+ 


,u>, 


n  ) 


dv' 


(14) 


where  in  two  and  three  dimensions,  the  imaginary  part  of  the  free  reference  space 
Green's  functions  are 


<  X  |  X ' »  i*i  ,  o  )  *  i  J0(a)ro) 


(15) 


and 
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1 


rMxjxJu.,  ')  = 


respectively,  and  where  3^=  |  X— X  * 


Taking  the  Hilbert  transform  of  (14)  with  respect  to  the  reference  slowness  yields 


1  (  (1  ( X ,  u  ,  z  '  ) 


'■  -j  if 


-  d  ■'  .  (x*  dv' 


+  j*i|  G^jxLxy^g 


m  r  11,2 

L  cL'(x' , 


7  ’  ■’  I  do'  *  ( X  ’  ,  ^  Jv  ' 


By  (15)  and  (16),  with  the  aid  of  [2,3]  and  twice  repeated  application  of  |4,5| 


it  fo I  lows  t hat 


on 

J  =  ,y 


X  X'  ,  , 


j  d0. 


3^w(u.r)  (x jx ’  ,  •  ■.  ’■ ' 


in  two  and  three  dimensions,  but  not  in  one  dimension.  Where  in  two  and  three 
dimensions  is  the  real  part  of  the  free  reference  space  Green's  function 


ryxix* , ..,.o  -  :  v 


and 


i  .  .  cos  ( uiro ) 

G  X  X'.w.o)  =  - 7 - 

r  '  4nr 


respect  i ve ly. 


Thus,  wi  th  the  aid  of  (18)  and  (19),  (17)  becomes 


^  — —  do'  =  j"  sgn(u>r)  G^(X|x’,ui,o )  p(X',u)  dv' 


+  I  sgn(ur)  G  (x|x',<*),o) - uc  a-’  I  <fr(x',u>)  dv' 

Jv  ^  Lc2(X',U)  J 


Restricting  (22)  to  positive  non-zero  frequencies  w,  permits  its  rewriting  as  the 
principal  value  integral  (i.e.,  excluding  integration  over  X=X') 


CO 

~f  dq>  *  P  J*G^x|x',u,o)  p(x',u)  dv' 


+  p  I  G  (X |x ' 


,(D,0)  ^ 


C2(x' ,w) 


-  w2  o2  I  (J;  (x'  ,w )  dv ' 


Since  the  imaginary  oart  of  the  free  reference  space  Green's  function  is  not  singular 
at  X=X It  follows  from  the  addition  of  (14)  and  (23)  that  for  positive  non-zero 
frequencies 
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1 


+  i  6  ( X ,  uj ,  o  ) 


*  P 


G ( X  | X  '  ,..:,0  ) 


(X' 


+  p  j  G(x  ;x* ,  uj,  i) 

*  v 


(24  1 


For  the  inverse  scattering  inverse  source  case  of  a  known  constant  medium  velocity  and 

unknown  sources  p(i.e.,  wave  equation  (1)),  (24)  reduces,  after  chosing  the  reference 

1  _ 

slowness  "  ,  for  to 

0  (X ,  to,  O  ) 

- ; -  do  +  i0(X,u), 

C  0 

At  that  reference  slowness,  the  direct  Kirchhoff  integration  (5)  of  the  wave  equation 
(1  )  can  be  wr i tten  as 


-)  = 


J 


G(X  X’ 


■a),-) 


r(x’,uG  dv 


(25) 


I 


(x,u.)  =  P  I  G  (X  I X  ’  i'(x’ 


■)  .i\ 


.(X. 


(241 


Since  in  two  and  three  dimensions  the  Green's  function  singularity  is  weak  and 
removable,  and  the  Kirchhoff  surface  integral  represents  the  incident  field. 

Thus  combining  (25)  and  (26)  yields  the  solution 


f  ( X ,  (o ) 


+  i  0  ( X  ,  .o 


,*■  . 1 X  ,  ,  0  ) 


(27) 


For  the  inverse  scattering  inverse  source  case  of  known  zero  sources  and  unknown 

medium  wave  velocity  c(X,u>i  i.e.,  wave  equation  (3),  (24)  reduces,  after  choosing  the  ' 

reference  slowness  o---,  for  ui>0,  to  j 

Co  3 
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■3(X,u),o  1 


f 


do  +  ) 


—  J 


— 


( x '  ,o) 


(X' , .  . 


1  V  ' 


(28) 


A  digression  concerning  the  wave  equation  (3)  is  now  in  order.  This  wave  equation 
can  be  r  ewr i 1 1 en  as 


(29) 


where  the  sources  p0  are  reference  sources  relative  to  the  arbitrarily  chosen  reference 
wave  velocity  c0,  given  by  the  relative  constitutive  equation 


(  3fl  ‘ 


and  »he  potential  V0is  a  reference  potential  relative  to  the  arbitrarily  choosen 
reference  wave  velocity  cot  given  by 


V„  = 


C-'  (X,uO 


(31  ) 


It  thus  follows  that  (28)  can  be  written  for  .  -0  as 


9(X,u),o) 

i 

-  -0 
Co 


do  +  iS(x,w 


f  >  i'otX’,- 


TT 


1  dv’ 


(32) 


At  this  reference  slowness  the  direct  Kirchhoff  integration  (5)  of  the  wave 
equation  (29)  can  be  written  as 


since  in  two  and  three  dimensions  the  Green's  function  singularity  is  w«ak  and  ruTcvahle, 
and  the  Kirchhoff  surface  integral  represents  the  incident  field. 

Thus,  combining  (32)  and  (33)  yields  the  (same  as  (27))  solution 


>-  c 


The  sources,  potential,  and  medi  i*n  wave  velocities  can  he  obtained  from  knowledge  of 
the  fields  in  a  variety  of  straight-forward  manners. 

One  might  be  tempted  to  attempt  to  simplify  the  solutions  (27,34)  by  applying  and 
executing  analytically  the  Hilbert  trans form  wi th  respect  to  the  reference  slowness 
dircectly  on  the  surface  integral  definition  (7)  of  the  effectal  field,  thus  obtain."- 
this  definition  in  terms  of  the  principal  value  of  the  real  part  of  the  free  rcfero'ur 
space  Green's  function,  instead  of  the  imaginary  part  of  this  free  reference  space  Green's 
function.  The  flaw  with  such  an  attempt  is  that  on  the  surface  of  integration,  awav  fro- 
the  singular  point  of  the  real  part  of  the  free  reference  space  Green's  function,  the 
principal  value  and  the  complete  singular  real  part  of  the  free  reference  space  Green's 
function  are  indistinguishable  and  identical,  and  the  application  of  Green's  theorem  as 
per  (7-12)  will,  by  the  differentiability  and  continuity  requirement  imposed  by  Green's 
theorem,  fail  to  generate  the  principal  value  of  the  real  part  of  the  free  reference  snace 
Green's  function  in  the  volume  interior  to  the  surface  of  integration,  which  would  have 
yielded  the  desired  solution,  but  generate  the  full  singular  real  part  of  the  free  reference 
space  Green's  function  in  this  interior,  \tfiich  fails  to  yield  the  desired  solution,  and 
yield  (a  slightly  modified)  version  of  the  integral  equation  (12),  instead  of  its  solution. 
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